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We study by Molecular Dynamics simulation the slowing down of particle motion in a two- 
dimensional monatomic model: a Lennard-Jones liquid on the hyperbolic plane. The negative 
curvature of the embedding space frustrates the long-range extension of the local hexagonal 
order. As a result, the liquid avoids crystallization and forms a glass. We show that, as 
temperature decreases, the single particle motion displays the canonical features seen in real 
glassforming liquids: the emergence of a "plateau" at intermediate times in the mean square 
displacement and a decoupling between the local relaxation time and the (hyperbolic) diffusion 
constant. 



1. Introduction 

Glass formation results from the slowing down of dynamics as a liquid is cooled to 
low enough temperature. Progress in a theoretical description of the phenomenon 
has been hampered by the shortage of an adequate and simple enough liquid 
model, which could play a role similar, say, to that of the Edwards- Anderson 
model in spin-glass theory [1]. Motivated by earlier work based on the concept 
of "frustration" O [3lll] , viewed as an incompatibility between extension of the lo- 
cally preferred order in a liquid and tiling of the whole space, we have recently 
studied a monatomic two-dimensional glassforming model: the Lennard-Jones liq- 
uid embedded in the hyperbolic plane [5j. The negative curvature characterizing the 
latter induces a frustration in the tendency of the system to order in a hexagonal 
structure (as it would do in the "flat" Euclidean space) and forces in an irreducible 
density of topological defects [6j. As a result, long-range ordering is thwarted and 
the two-dimensional monatomic liquid can now form a glass. Elsewhere, we have 
shown that (1) the "fragility" of the liquid, i.e. the extent to which the temper- 
ature dependence of the relaxation time deviates from a simple Arrhenius law 
and displays "super-Arrhenius" behavior, can be tuned at will by changing the 
curvature of the embedding space and (2) that the dynamics become spatially het- 
erogeneous as temperature decreases f5^. A potential interest of the model lies in 
its two-dimensional and monatomic character: this opens the way to theoretical 
developments based on the presence of well-identified point-like topological defects 

mm- 

In this short article, we consider another aspect of the dynamics of the Lennard- 
Jones liquid on the hyperbolic plane. We monitor the single particle (translational) 



'Corresponding author. Email: sausset@lptmc.jussieu.fr 

ISSN: 1478-6435 print/ISSN 1478-6443 online 

© 200x Taylor & Francis 

DOI: 10.1080/1478643YYXXXXXXXX 

http:/ /www. informaworld.com 



December 19, 2008 



11:37 Philosophical Magazine Andalo 



2 



motion as one lowers the temperature and approaches glass formation. More specif- 
ically, we study the properties of the long-time "diffusive" behavior and the decou- 
pling of diffusion from local structural relaxation that one often associates with a 
breakdown of the Stokes- Einstein r elation [5^, 'n]. One motivation for this study is to 
verify that the hyperbolic metric does not introduce spurious effects and that the 
present model indeed displays the phenomenology observed in real glassforming 
liquids. 

2. Model 

We have considered what seems to be the simplest microscopic liquid model 
that is subject to "geometric frustration": a monatomic liquid on the hyper- 
bolic plane. Atoms interact through the standard Lennard-Jones pair potential 
v{r) = 4e ((<7/r)^^ — (u/r)^), where the distance r is defined with the hyperbolic 
metric. The hyperbolic plane (denoted H^) is a Riemannian surface of constant 
negative curvature {K < 0), whose metric in polar coordinates (r, (p) is given by 



where k is related to the Gaussian curvature K of the plane by K = —k^. 

The Euclidean counterpart of this system (k = 0) orders to a hexagonal "crys- 
tal" via either a weakly first-order transition or a sequence of two nearby contin- 
uous transitions separated by an intermediate "hexatic phase" [TJ |T0] . However, in 
nonzero curvature («; 7^ 0), the hexatic/hexagonal order that occurs in flat space 
can no longer tile the whole plane. If the curvature, as characterized by the dimen- 
sionless parameter ko", is small enough, the hexagonal order is still locally favored. 
Then, the local and global orders are incompatible and the system is (geometri- 
cally) frustrated O El Eli therefore plays the role of an intrinsic frustration 
length. The ordering transition that takes place in zero curvature at (or around) 
a temperature T* is avoided [11 for any nonzero curvature, no matter how small; 
the system remains in the liquid phase even below T* and may subsequently form 
a glass[5]. By an abuse of language, we shall refer to the liquid below T* as be- 
ing "supercooled" even though it does not compete with a crystalline phase when 
the curvature is nonzero. (As discussed in Refs. Ul [11], one can associate T* in 
real three-dimensionnal glassforming liquids with a temperature that is empirically 
found to mark the onset of "anomalous" behavior characterized by super- Arrhenius 
temperature dependence of the relaxation time and non-exponential time depen- 
dance of the relaxation functions [12j.) 

We have performed a Molecular Dynamics simulation of the model. Details on the 
procedure that requires nontrivial extensions of the simulation method commonly 
used in Euclidean space are given in Ref. |5j. The simulation is performed in the 
microcanonical ensemble (A^, V, E) and octagonal periodic boundary conditions are 
used [3 [13]. The control parameters are the temperature Te~^, the density pcr^, and 
the frustration associated with space curvature na. 

3. Self-diffusion on the hyperbolic plane 

In this article, we focus on several aspects of the single particle motion. In a real 
liquid, it is well known that the long-time limit of such a motion is described by a 
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diffusion process cliaracterized by tlie (self-) diffusion coefficient D. Standard calcu- 
lations show that the mean square displacement of the particles goes asymptotically 
as predicted by the Einstein relation: (|r(t)^^ ~ ADt, where r{t) = |r(t) — r(0)| is 
the absolute value of the displacement of a given particle during the time t. In a 
curved space, this relation does not hold anymore and one may wonder what is 
the generalization of the diffusive limit. To answer this point, we first take a little 
mathematical detour for considering the Brownian motion in the hyperbolic plane. 

The Brownian motion in the hyperbolic plane has recently been studiedjl4| 115] . 
It is based on a Fokker-Planck equation describing the time evolution of the prob- 
ability density P{r, t) associated with finding the particle at a distance r of its 
starting point after a time t: 



dP{r, t) 
dt 



D- 



d 



dP{r, t) 
dr 



K coth{nr)P{r, t) 



(2) 



whose solution is 



P{r,t) = — -= 3- sinh(Kr) / dy , (3) 

2\/2Tr{Dt)^ Jr -\/cosh(Ky) — cosh(Kr) 

where D is the (hyperbolic) diffusion coefficient. Equation ([3| allows one to calcu- 
late the time dependence of the mean square displacement^r^(t)) and the mean 
(absolute) displacement {r{t)). An asymptotic analysis then shows that when the 
distance travelled r ^> ballistic diffusion is obtained with (|r-^(t)) ~ {Dt)"^ 

and (r(t)) ~ Dt^ whereas when r <^ k^^, ordinary (Euclidean-like) diffusion is 
found with (r^(t)) ~ Dt and {r{t)) ~ (Dt)^/'^. (Note that, of course, it is the same 
diffusion coefficient D that appears in these two regimes.) 

The working assumption we make to generalize the treatment of liquids in ordi- 
nary Euclidean space is that the single particle motion in the liquid embedded in 
the hyperbolic plane converges at long times to a process described by the above 
Brownian motion. 

In practice, and to increase the domain over which one can use a simple fitting 
formula to extract the self-diffusion coefficient D from the single particle motion, 
it is convenient to devise a function of the (absolute) displacement, /(r(t)), that 
behaves similarly in the intermediate Euclidean-like diffusive regime and in the 
asymptotic ballistic diffusive regime, i.e. such that {f{r)) ^ Dt over a wide time 
domain. An expression of / can be obtained by using the Fokker-Planck equation 
([2]) and it reads 



^/ / 9, /I -|- cosh(Kr(t)) \ 
f{r{t))=K-Hxi[^ J • 

One easily checks that indeed {f{r)) = Dt for a Brownian particle. Thus, by 
calculating {f{r{t))) = (1/A^) Yld=i /('"iC*)); where the sum is over all the atoms of 
the system, and comparing it with a linear time dependence, one can verify that the 
atomic motion in the Lennard-Jones liquid becomes diffusive at long enough times 
and obtain the coefficient of self-diffusion D without worrying about the crossover 
between ordinary and ballistic diffusion. That it is actually the case is illustrated 
on the model with frustration parameter na = 0.2, at the density pa^ = 0.91, 
and at a fairly high temperature (T = 2.17, much above the "avoided transition" 
temperature T* ~ 1.3) so that the atoms move distances much larger than k""^ in 
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Figure 1. Log-log plot of the mean absolute displacement in units of and of {f{r{t))) (see Eq. 

|j4]|). The dashed line has a slope equal to 1 and the dotted-dashed one a slope of 1/2. {f(r{t))) has a 
roughly linear time dependence at all times, which corresponds to a diffusive regime in the hyperbolic 
sense. The parameters are pa^ = 0.91, T = 2.17, and kct = 0.2. The system with octagonal periodic 
boundary conditions comprises 287 atoms. 

the simulation time. The results can be seen in figure [T] where the mean absolute 
displacement {r{t)) and the associated {f{r{t))) are plotted versus time. One can 
see that, as expected, exhibits two different regimes in the diffusion limit: 

one following ordinary diffusion for r < and one following ballistic diffusion for 
r S> thus sets the scale at which the crossover between the two regimes 

takes place. The figure also illustrates that it is easier to extract the self-diffusion 
coefficient D from fitting (/(r(t))), which behaves as Dt over a very wide domain. 

4. Single particle motion in the supercooled liquid 

In the "supercooled" liquid (i.e. the liquid below T* , see above), one observes a 
dramatic slowing down of the dynamics. As shown in Ref. [5J, the super Arrhenius 
character and the associated fragility increase as one decreases the frustration Ka. 
Here, we focus on a moderately small frustration, Ka = 0.1, and on two features 
of the single particle dynamics at sufficiently low temperatures: the emergence of 
a "plateau" in the mean (absolute) displacement of the atoms and the decoupling 
of the diffusion and the local structural relaxation. 

4.1. Plateau in the mean (absolute) displacement 

We now compute the mean absolute displacement (r(t)) of the atoms in a system of 
density pa'^ = 0.91 for the liquid phase both above and below the avoided transition 
temperature T* ~ 1.3 at which ordering occurs in the Euclidean plane: see Fig. 
[2} At high temperature above T*, (r(t)) is very close to its counterpart in the 
Euclidean plane in the domain covered, i.e. for (r) up to a few atomic diameters. 
On such distances that are significantly smaller than the intrinsic frustration length 
K"^ (recall that k,~^ = Wa here), the dynamics is dominated by local relaxation 
processes unaffected by the curvature. At short times, {r{t)) is ballistic due to 
the Newtonian dynamics of the atoms and it crosses over to a diffusive behavior at 
longer times. When decreasing the temperature and entering the supercooled liquid 
regime, a plateau appears in (r(t)) at intermediate times and for a displacement of 
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Figure 2. Log-log plot of the mean absolute displacement {r{t)) of the atoms in units of the atomic 
diameter cr. For this case, Ktr = 0.1, which means that when (r) ^ tr, it is only a tenth of . Three 
temperatures are represented: from top to bottom, T = 1.53, T = 0.96, and T = 0.47, with the avoided 
transition temperature being at T* ~ 1.3. As T decreases, a "plateau" emerges at intermediate times and 
it is well developed at the lowest temperature. At longer times, a diffusive behavior is found, as shown by 

the dashed line whose slope is equal to 1/2. 

a fraction of a. This plateau becomes more pronounced as T is lowered and for the 
lowest temperature {r{t)) stays roughly constant for about two decades of time. 
This phenomenon is usually described ge effect": particles are trapped by 

their neighbors (which are themselves trapped) and they can only rattle before 
being able to escape the cage and change their neighborhoodjT6]. At long times, 
the atomic motion becomes diffusive. All of this is very similar to what is found in 
simulations and experiments on 3-dimensional glassforming liquids|T71 118j. 

4.2. Decoupling of diffusion and local structural relaxation 

Another intriguing experimental observation found in supercooled liquids is the 
decoupling of the characteristic time scales measured by different probes (transla- 
tional diffusion on one hand, viscosity and a-relaxation time on the other hand) and 
the associated breakdown of the Stokes- Einstein relation [8l [9]. In ordinary liquids 
(above the melting point), a homogeneous description is found to apply down to 
molecular scale and the hydrodynamic-based Stokes-Einstein relation, ^ = const. 
with D the self-diffusion coefficient and rj the viscosity, is well obeyed |19]. However, 
experimental measurements have shown that this relation does not hold anymore 
in the supercooled liquid regime, where the ratio ^ increases with decreasing 
temperature by up to several orders of magnitude [H |9l [20] . On the other hand, the 
a-relaxation time Ta obtained from rather local measurements, i.e. probes involving 
particle displacements over a small, molecular scale, such as rotational and struc- 
tural relaxation measurements, follow rather well the temperature dependence of 
the viscosity and therefore also decouple from the self-diffusion process that implies 
particle translational motion over large, super-molecular distances. 

To study such a potential decoupling in our model, we have monitored both the 
self-diffusion coefficient D (see above) and the a-relaxation time obtained from 
the hyperbolic generalization of the self-intermediate scattering function Fs{k,t), 
when studied for a wavevector k ~ close to the first maximum of the static 
structure factor [5i. The temperature dependences of T/D and Tq are displayed in 
the inset of Fig.|3]in an Arrhenius plot. As can be seen from the plot, the deviation 
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Figure 3. Decoupling ratio between the diffusion and the alpha-relaxation versus the inverse 

temperature. The ratio is normalized to 1 at high temperature and the temperature is scaled by T*. In 
the inset we separately display and Tc in an Arrhenius plot; the scales on the j/-axis are shifted in 
order to bring the two sets of data together at high temperature. 



from Arrhenius behavior is not very marked for this cm'vature (kc = 0.1) and 
this density {pa'^ = 0.91); the corresponding glassforming hquid is rather "strong" 
(or in a terminology that is somewhat more appropriate for atoms interacting via 
weak nondirectional pair potentials, has a small fragility [2T]). Nonetheless, there 
is a significant increase, by an order of magnitude, of the analog of the Stokes- 
Einstein ratio, temperature decreases in the supercooled regime: see 

Fig. |3j One expects that the magnitude of the decoupling would increase as the 
fragility of the liquid increases. This can be investigated in our model by decreasing 
the curvature kcj of the embedding space, and this work is in progress. 



5. Conclusion 

We have illustrated that the frustration-based microscopic model of a two- 
dimensional monatomic Lennard-Jones liquid embedded in the hyperbolic plane 
displays many features of the phenomenology of supercooled liquids (see also [5] ) . 
In this short article, we have focused on the single particle motion, showing that 
the long-time behavior is described by a hyperbolic generalization of the diffu- 
sion process and devising an efficient way to extract the self-diffusion coefficient. 
The two canonical features of glassforming liquids as temperature is lowered that 
we have retrieved are the emergence of a plateau in the time dependence of the 
mean (absolute) displacement of the atoms and a decoupling of the diffusion and 
the a-relaxation. One main advantage of the model is that by changing the space 
curvature one varies the fragility of the system so that one can envisage the phe- 
nomenology of glassforming liquids at constant molecular parameters but varying 
fragility. In addition, one expects that the model may lend itself to a theoretical 
treatment based on the consideration of the point-like topological defects that are 
characteristic of two-dimensional, one-component systems. 
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